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ABSTRACT 



We deal with the effects induced on the orbit of a test particle revolving 
around a central body by putative spatial variations of dimensionless fundamen- 
tal coupling constants £. In particular, we assume a dipole gradient for C( r )/C 
along a generic direction k in space. We analytically work out the long-term vari- 
ations of all the six standard osculating Keplerian orbital elements parameterizing 
the orbit of a test particle in a gravitationally bound two-body system. Apart 
from the semi-major axis a, the eccentricity e, the inclination /, the longitude of 
the ascending node Q, the longitude of pericenter w and the mean anomaly Ai 
undergo non-zero long-term changes. By using the usual decomposition along 
the radial (R), transverse (T) and normal (N) directions, we also analytically 
work out the long-term changes AR, AT, AiV and Avr, Avt, Avn experienced 
by the position and the velocity vectors r and v of the test particle. Apart from 
AiV, all the other five shifts do not vanish over one full orbital revolution. In the 
calculation we do not use a-priori simplifying assumptions concerning e and I. 
Thus, our results are valid for a generic orbital geometry; moreover, they hold for 
any gradient direction k. We compare our predictions to the latest observational 
results for some of the major bodies of the solar system. The largest predicted 
planetary perihelion precessions occur for the rocky planets, amounting to some 
10 -2 — 10~ 3 milliarcseconds per century. Apart from the Earth, they are 2 — 3 
orders of magnitude smaller than the present-day accuracy in empirically de- 
termining the corrections Am to the standard Newtonian-Einsteinian planetary 
perihelion rates. Numerically integrated time series of the interplanetary range 
for some Earth-planet pairs yield Stark-like signatures as large as 0.1 — 10 mil- 
limeters; future planned planetary laser ranging facilities should be accurate at 
a cm level. The long-term variations of the lunar eccentricity and perigee are of 
the order of 10 -14 yr -1 and 10 -4 milliarcseconds per year, respectively, while the 
change AR in the radial component of the Moon's geocentric orbit is as large 
as 0.8 microns per orbit. A numerically calculated geocentric lunar range time 
series has a maximum nominal peak-to-peak amplitude of just a few millimeters, 
with an average of 0.3 microns over 30 yr. The present-day accuracies in deter- 
mining e and w for the Moon are 10 -12 yr -1 and 10 _1 milliarcseconds per year, 
respectively. The APOLLO facility should be able to determine on a continuous 
basis the Earth-Moon range with a millimeter accuracy. 

Subject headings: gravitation; celestial mechanics; ephemerides; planets and satellites: 
general; Moon 
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Introduction 



Testing constancy of the fundamental constants entering the basic laws of physics is of 
the utmost importance for our understanding of the nature of the gravitational interaction 
and of the doma i n of v alidity of the Einsteinian general theory of relativity; for a recent 
review see 



Uzanl ( 1201 ll ) and references therein. 



Since the matter-energy content U = mc 2 of material bodies generally depends on 
the parameters of the Standard Model, a spatial variation in one of them will induce an 
extra-force on a body of mass m 



-W 



where3 £ is an adimensional fundamental parameter like, e.g., the fine structure constant or 
the electron-to-proton mass r atio, and c is the speed of li ght in vacuum. In particular, for a 
dipole-type spatial variation (IDamour fc Donoghuell201ll ) 



c 



1 + B k-r 



of £, the force is (IDamour fc Donoghuell201ll ) 

F = -mQBc 2 k, 

in which 



Q 



C dm 



(2) 

(3) 
(4) 



m d( 

is a dimensionless "charge". In eq. (j2J) and eq. (J3J) B is a slope parameter, having 
dimensions of L _1 , relative to a direction in the space determined by the unit vector k. For 
example, for the same direction^] 



k = {0.50 ± 0.08, -0.19 ± 0.04, -0.84 ± 0.04} 



with respect to an ecliptic frame, it was found ( jWebb et al.ll2010l ) 

B = (1.10 ± 0.25) x 10^ 6 Glyr^ 1 = (1.16 ± 0.26) x 10~ 31 m" 1 



(5) 



(6) 



lr The dimensionless ratios o f various fundamental parameters are dubbed r iy i = 1,2,... 
by lDamour fc Donoghud (120 111 ). 

2 It correspond s to equatorial coordinates RA= 17.4 ± 0.6 hr, DEC= —58 ± 6 deg 
flWebb et alJl2010h . 
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(7) 



for the fine structure constant, and (IBerengut et al.ll201ll ) 

B = (2.6 ± 1.3) x 10~ 6 Glyr" 1 = (2.7 ± 1.3) x 10~ 31 m" 1 

for the electron-to-proton mass ratio. 

If Q is not the same for all bodies, then a non-zero, net relative acceleration of 
Stark-type 

A = A B -A A = -AQBc 2 k, (8) 

where 

AQ = Q B - Q A , (9) 

occurs for a two-body system A-B. Notice that eq. (jHJ), which implies a violation of the 
equivalence principle, holds for a generic adimensional parameter (; in principle, the total 
extra-acceleration is the sum of all the terms like eq. (JSJ) due to the gradients of the various 
(. As far as the magnitude of eq. (|S]) is concerned, it is certainly quite small. Indeed, from, 
say, eq. (jBJ) it is 

A < 10" 14 m s^ 2 ; (10) 

suffices it to say that the standard Newtonian inverse-square law for the Sun-Earth system 
yields 

^Newton - ID" 3 m S 2 . (11) 

The plan of the paper is as follows. In Section [5] we analytically work out the long-term 
effects caused by a Stark-type extra-acceleration of the form of eq. (JSJ) on the motion 
of a test particle orbiting a central body. We do not make any a-priori assumption 
concerning both the particle's orbital geometry and the direction k. Consequences of a 
vio lation of the equiva l ence p rinciple referred to a fixed direction in space were investigated 
by ID amour fc S chafer! (119911) in the framework of binary pulsars, and preliminarily by 
Damour fc Donoghuel ( 2011 ) for the Earth-Moon system. In Section [3] we compare our 
results to the latest empirical determinations from solar system observations. Section H] is 
devoted to the summarizing our results. 



2. Analytical calculation of the orbital effects 

The standard Keplerian orbital elements of the orbit of a test particle are the 
semi-major axis a, the eccentricity e, the inclination /, the longitude of the ascending node 
Q, the argument of pericenter u, and the mean anomaly M.. While a and e determine the 
size and the shaped, respectively, of the Keplerian ellipse, /, f2, cu fix its spatial orientation. 



3 The eccentricity e is a numerical parameter for which < e < 1 holds; e = corresponds 
to a circle. 
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/ is the inclination of the orbital plane to the reference {x, y} plane, while Q is an angle 
in the {x, y} plane counted from a reference x direction to the line of the nodes, which is 
the intersection of the orbital plane with the {x, y} plane. The angle co lies in the orbital 
plane: it is counted from the line of the nodes to the pericenter, which is the point of closest 
approach of the test particle to the primary. In planetary data reduction the longitude of 
the pericenter w = Q + u is customarily used by astronomers, although it is a "dogleg" 
angle. The argument of latitude u = u + / is an angle in the orbital plane which reckons 
the instantaneous position of the test particle along its orbit with respect to the line of the 
nodes: / is the time-dependent true anomaly. The mean anomaly is defined as 



M=n{t- t 



VII 



'121 



where 



n 



\jGM/a? 



(13) 



is the Keplerian mean motion related to the Keplerian orbital period by n = 2ir/P^, and t p 
is the time of passage at the pericenter. In the unperturbed two-body pointlike case, the 
Keplerian ellipse, characterized by 



x = r (cos Q cos u — sin Vt sin u cos I) 
^ y = r (sin Q cos u + cos Q sin u cos /) 
z = r (sin u sin /) , 



(14) 



and 



a{\ -e 2 ) 



(15) 



1 + e cos / 

neither varies its shape nor its size; its orientation is fixed in space as well. 

A smalfl perturbing acceleration A, like eq. (JSJ), of the dominant inverse- square 
Newtonian term Aj<( ewton induces slow temporal changes of the osculating Keplerian orbital 
elements. The Gauss equations for their variation, valid for any kind of acceleration 



4 Cfr. eq. (ITU]) and eq. f lTTj) . and the figures for AQ in Table [TJ below. 
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whatever its physical origin may be, are ( jBertotti et al.ll2003l ) 

[eA R smf + A T (*;)], 



f da 

dt 



dt 



dl 
dt 



dQ 
dt 



dt 



^{A R smf + A T [coB/+Hl-S)]} 

7= bAjv ( L ) COSW, 



1 



? A N (I) sin 



(16) 



a. 



na sin 7%/l— « a 

-Aflcosj + Ar (l + £)sin/ 



+ 2 sin 2 (|) dQ 



dt ' 



n-^(S)-vT^(f + cos/f). 

In eq. f|T6|) p = a(l — e 2 ) is the semi-latus rectum, and Ar, At, A^ are the radial, transverse 
and out-of-plane components of the disturbing acceleration A, respectively. In a typical 
first-order perturbative calculation, they have to be computed onto the unperturbed 
Keplerian ellips^E according to 



;i7) 





= A 


R, 


< A T 


= A 


T, 




= A 


•N, 



where the unit vectors along the radial, transverse and out-of-plane directions are 

cos f2 cos u — cos / sin Vl sin u, 

R = < sinficosu + cos/cosf2sinu, 

sin / sin u, 
— cos f2 sin u — cos / sin f2 cos u, 

T={ — sinfisinw + cos / cos Q cos it, 

sin / cos u, 



(19) 



5 In principle, one may choose a different reference orbit including general relativity itself 
( iCalura et al.lll997lll998l ). However, completely negligible mixed Einstein-Stark terms would 
result. We will not deal with them. Such a conclusion will be a-posteriori confirmed by the 
extremely small magnitudes of the presently computed Newton-Stark effects (See Section [3] 
below). 
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sin/sinf2, 
N = < — sin/ cos il, 
cos/. 

V 

The result is, after some algebra, 

Ar = Bc 2 AQ < — cosu ( k x cosf2 + k y smil ) — sinw 



(20) 



k z sin / + cos / ( k v cos — k x sin £7 



^ A T = Bc 2 AQ 



—k z cos u sin / — cos / cos u ( /c„ cos Vt — k x sin Vt ) + sin w ( k x cos f2 + k\ sin O 



A' 



-Bc 2 AQ 



k z cos / + sin / ( k x sin Vt — k v cos 



(21) 



In the case of eq. ([H]), it turns out that it is computationally more convenient to use 
the eccentric anomaly E instead of the true anomaly /. Basically, E can be regarded as a 
parametrization of the usual polar angle 9 in the orbital plane, being defined as 

M = E - esinE. (22) 

To this aim, useful conversion r elations, to be used when eq. (I2TT) is inserted in the 
right-hand-sides of eq. ( 1TB|) . are (IBertotti et al.ll2003l ) 



' cos ,/ 

sin/ 
r 

dt 



cos E—e 
1— ecos_B' 



e 2 sini? 
1— ecos_B ' 



a(l — ecosE), 



(23) 



Cumbersome calculation yield the variations of all the Keplerian osculating orbital elements 
averaged over one orbital period as 



~dt 



n 

27 



— — dt, * = a, e, /, Q, w, M, 
dtJ G 



(24) 



where it is intended that eq. (I2"3"j) is used in the integrand of eq. ( 12~4"|) and itself: (d^/dt) G 
represents the right-hand-sides of eq. ffl6|) for the generic orbital element Note that, as 
usual in perturbation theory, a, e, /, Q, u are kept fixed in performing the integration with 
respect to the fast variable E. Indeed, their actual variations due to several non-Keplerian 
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effects occur over characteristic timescales which are quite longer than P^. In the following, 
the brackets (• • • ) denoting the average over one orbital period will be omitted for brevity. 
One finally gets 



f da 

dt 



de 
dl 



dJ_ 

dt 



dn 

dt 



IF 



0. 



3Bc 2 AQV~r^P 
2an 



k z sin / cos u + cos I cos to ( k y cos Vt — k x sin Q 



— sin cu (k x cos Q + k y sin Q j 



3Bc 2 AQecosu 
2an\J\~ e 2 



3Bc 2 AQesinu 
2an\J\— e 2 



k z cos 7 + sin I ( k x sin VL — k v cos f2 



^ sin Q — k v cos f2 + /c z cot / 



(25) 



2^0? { ( e2 " !) c °s w cos + k y sin 



+ sin a; 



— k z sin / + (e 2 — cos 7) cos — /c x sin j + e 2 k z tan (|) | , 



dM 
dt 



3Sc 2 AQ(l+e 2 
2ean 



cos cu i k x cos VL + k y sin f2 j + 



+ sinw 

v 

We remark that the expressions in eq. ([25]) are exact in the sense that no simplifying 
approximations either in e or in I were assumed in the calculation; moreover, they are 
valid for a generic direction k of the dipolar gradient of (. It can be noticed that the 
semi-major axis remains unchanged, while the long-term variations of the inclination and 
the node vanish for circular orbits. The formula for dQ/dt becomes singular for / — > 0; the 
same occurs for dw/dt and dAi/dt as well for e — > 0. In general, the long-term changes 
of eq. fl25|) are not secular trends because of the modulations introduced by the slowly 
time-varying orbital elements themselves occurring in real astronomical scenarios like the 
Earth and the Moon, and the Sun and its planets. In the calculation yielding eq. ( 1251) it 
was assumed that their frequencies were much smaller than the orbital one, so to keep them 
constant over one orbital revolution. 

The instantaneous changes of the R — T — N components of the test particle's position 



k z sin I + cos / (k y cos Q — k x sin Qj 
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vector r can be worked out from the following general expression (jCasottolll993l ) 
' AR(f) = (9Ao(/)-ocoB/Ae(/)+(^)AA<(/) I 

< AT(f) = asin/Jl + ^y] Ae(f) + r[Au(f) + cos IAQ(f)} + ^ VT^AM(f), 
AN(f) = r[smuAI(f) - cosm sin/Afi(/)] , 

(26) 

In the case of eq. (jHJ), we have that the R — T — N shifts of the position, averaged over one 
orbital revolution, are 



AR 



3irBc 2 AQ\/l-e 2 

r,2 



k z cos oj sin I + cos I cos u ( k v cos Vt — k x sin f2 



S AT 



sin u ( fc x cos Vt + k y sin f2 



(27) 



+ sinw 



AN = 0. 



/c* sin I + cos / ( k y cos f2 — k x sin f2 



Also the expressions of eq. ( 127|) are exact in both e and /; notice also that they present 
no singularities for both e — > and I — > 0. Moreover, they, in general, vanish neither for 
circular orbits nor for / = 0. 

For the instantane ous R — T — N perturbations of the velocity vector v we have, in 
general, (jCasottolll993[ ) 



Av R (f) 

Av T (f) 
Av N (f) 



n sin / 
Vl-e 2 



|Aa(/) + ^Ae(f) 



[Aw(/) + cos JAO(/)] - (s£ AM(/) 



nav l— e 



, Aa(/) + ^5p#Ae(/) + ^ [A W (/) + cos JAfl(/)] 
[(cos u + e cos w) AI(f) + (sin u + e sin a;) sin 7Afi(/)] . 



(28) 

In the case of eq. (J5]), eq. ( 128]) yields the following expressions for the velocity changes over 
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one orbital revolution 



< Av T 



Av N 



3itBc AQ[l+e(2 e)] ^ cQg ^ ^ ^ Q _|_ ^ g j n ^ 



+ sinw 



/c z sin I + cos / (k y cos f2 — &; x sin j 



3ttBc 2 AQ 
(l-e)n 



fc 2 cos cu sin J + cos / cos ui ( fc v cos VL — k x sin fl 



(29) 



sin a; cos fl + k y sin fl 



'iivBc 2 AQe 
(l-e)n 



fc z cos / + sin I ( k x sin Vt — k y cos 



Also the expressions of eq. ( 1291) are exact in both e and J, and are not singular for any 
particular value of them. Notice that Avn vanishes for circular orbits. 



3. Confrontation with the observations 

For the sake of definiteness, in the following we will focus on the fine structure constant 
whose spatial variation has the most stringent empirical support so far. 

The following general considerations are in order to properly contextualize and 
understand the content of the next sections. Our aim is to evaluate as realistically as 
possible the magnitude of the predicted anomalous effects on quantities which are empirically 
determined from observations by astronomers. It is important to note that they were 
produced by adopting force models including all the standard known Newtonian/Einsteinian 
dynamics which, thus, did not include the dipolar effects considered here. In principle, 
such observationally determined quantities should account for them, but it is quite possible 
that in the estimation procedure of, say, the planetary state vectors the exotic signatures 
were partially or totally absorbed, especially if their magnitude is small with respect to 
the accuracy of the orbit determination process. Thus, the subsequent comparison of the 
theoretically computed Stark-like effects with the present-day empirical accuracies is just 
a preliminary, although quite reasonable and necessary, step to investigate if the putative 
orbital anomalies considered here may be considered as potentially measurable in further, 
dedicated analyses. We do not intent to put on the test their actual measurability: indeed, 
to this aim it would be necessary to reprocess the entire lunar/planetary data sets with 
ad- hoc modified dynamical models explicitly including the Stark-type effect of interest, and 
dedicated solve-for parameter(s) should be estimated in the new global fit. Of course, such 
a non-trivial task would be worthwhile only if there are some hopes that the effects one is 
interested in are, at least indicatively, detectable. 
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3.1. The planets 



In order to exploit the latest results in the field of planetary orbit determination, 
in Table Q] we quote the value s of AQ for several Sun-plan ets pairs. We computed 
them by means of eg. (42) in [5 amour fc Donoghud (120 111 ), based on earlier results 
( Damour fc Donoghue 2010al Jbl). As it can be expected, AQ is larger for the rocky 
planets, mostly made of iron cores and silicate mantels, with respect to the gaseous giants 
l ike Saturn, for wh i ch accurate orbital data ar e now available from the Cassini mission 
(IFienga et al.ll2010l iFiengal |2010| ; IPitjeval 120101 ) . whose composition (IFouchet et al.1 120091 ) 
is more similar to that of the Sun. For our parent star we ass umed that hydrogen and 
helium constitute 74% and 25% of its total mass, respectively (jGrevesse &: Sauvallll998l ). 
Concerning Mercury, it is known since a long time (jUreylll95ll . Il952l ) that its iron-to-silicate 
ratio A must be much larger than that of any other terrestrial planet and the Moon; for 
it we assumed A = 70 : 30 (IBenz et al.ll2007l ). Venus has a rather similar composition 
with respect to the Earth (|Lewisll2004l ). for whic h we a s sume A = 32 : 68 (jLewisI 120041 ) . 
According to a spacecraft data-based analysis by lYoderl (119971 ) . the venusian iron/silicate 
ratio ifl A = 26 : 74. Mars is an intermediate case between the Earth and the Moon 



(jLewisI 120041 ). which is almost entirel y made of silicates (jWilliams &: Boggsl 120091 ); latest 
determinations from orbiting probes ( Konopliv et al .1 1201 lh allow to infeit] A = 21 : 79 for 
the red planet. For Saturn we take a mass fraction for h elium of the order of 0.13 from the 
latest Cassini-based measurements (IFouchet et al.l 120091 ). while the rest is assumed to be 
made mostly of hydrogen. 



Table 1: Values AQ = Q p — Qq for various planet-Sun pairs inferred from eq. (36) of 
Damour fc Donoghud (1201 ll ) for the case of the fine structure constant. See the discussion 



in the text for the details concerning the composition of the various bodies. 



Mercury 


Venus 


Earth 


Mars 


Saturn 


1.979 x 10~ 3 


1.539 x 10~ 3 


1.599 x IO- 3 


1.489 x lO" 3 


-2.8 x 10- 5 



6 Indeed, according to Figure 2 by lYoderl (119971 ). the mean core density of Venus is g c 
10.16 g cm -3 for a value of the core radius of R c = 3100 km. 

7 Cfr. with Figure 26 of (jKonopliv et al.l 120111). We used p r 
density and R c 



6.7 g cm 

1680 km for the core radius (jKhan fc Connolly! 120081 ) . 



for the core 



3.1.1. The perihelion precessions 



By using the figures of Table [T] and the result of eq. fl25|) , assumed to be valid in a 
heliocentric reference frame with the mean ecliptic and equinox at the epoch J2000.0, we 
compute the long-term perihelion precessions induced by eq. ([8]). We adopt eq. ([5]) for k 
and eq. ([6]) for B. Table [2] displays our results in milliarcseconds per century (mas cty -1 ). 

Table 2: Perihelion precessions w, in milliarcseconds per century (mas cty -1 ), for the four 
inner planets and Saturn induced by a dipole-type gradient of the fine structure constant 
according to eq. fl25|) . We used eq. ([5]) for k, eq. for B, and the figures of Tabled] 
for AQ. The Keplerian orbital elements of all the planets were kept fixed to their values 
referred to the mean ecliptic and equinox at the epoch J2000.0: they were retrieved from 
the NASA WEB interface HORIZONS (Author: J. Giorgini. Site Manager: D. K. Yeomans. 
Webmaster: A. B. Chamberlin). 



Mercury 


Venus 


Earth 


Mars 


Saturn 


-0.0016 


-0.0369 


-0.0244 


0.0017 


0.0004 



They should be compared with the latest empirical determinations of the corrections 
Arij to the standard planetary perihelion precessions obtained by fitting accurate dynamical 
force models, which include most of the standard Newtonian and Einsteinian dynamical 
effects, to observational data records spanning about one century. They are resumed in 
Table EJ 

It can be noticed that the predicted precessions of Table [2] are, in general, smaller than 
the present-day uncer tainties in est i matin g Aw by 2 — 3 orders of magnitude, apart from 
the value obtained by iFienga et al.l (120101 ) for the Earth's perihelion rate by including some 
VLBI points for Jupiter. 



3.1.2. The Earth-planet ranges 

Since a direct, unambiguous observable is the range p between the Earth and a planet^, 
and in view of the po ssible future implementation of interplanetary laser ranging facil ities 
accurate to cm level (ITuryshev fc Williams! 120071 : iMerkowitz et al.l 120071 : iDegnanl 120081 ) , we 
investigate the impact that a violation of the equivalence principle due to a Stark-type 



8 It is particularly true when a spacecraft specially outfitted for ranging measurements 
orbits the planetary target of interest. 
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Table 3: Estimated corrections Am, in milliarcseconds per century (mas cty 1 ), to the 
standard Newtonian-Einsteinian secular precessions of the longit udes of the p erihelia m 
of the e ight planets plus Pluto determined wi th the EPM2 008 f lPitieval lioioh . the IN- 
POP08 flFienga et alJlboioj ). and the INPOPlOa flFiengallioioh ephemerides. Only the usual 
Newtonian-Einsteinian dynamics was modelled, so that, in principle, the corrections Am 
may account for any other unmodelled/mismodelled dynamical effect. Concerning the val- 
ues quoted in the third c o lumn from the left, they correspond to the smallest uncertainties 
reported by lFienga et al.l (120101 ) . Note the small uncertainty in the correction to the preces- 



sion of the terrestrial perihelion, obtained by processing Jupiter VLBI data (IFienga et al. 
2010h . 



Planet 


Am CPitieva 2010) 


Am fFienea et al. 2010) 


Am fFienga 2010) 


Mercury 


-4 ±5 


-10 ±30 


0.2 ±3 


Venus 


24 ±33 


-4 ±6 




Earth 


6±7 


0± 0.016 




Mars 


-7 ±7 


0±0.2 




Jupiter 


67 ±93 


142 ± 156 




Saturn 


-10 ± 15 


-10 ±8 


± 2 


Uranus 


-3890 ± 3900 


± 20000 




Neptune 


-4440 ± 5400 


± 20000 




Pluto 


2840 ±4510 







acceleration like eq. flS} may have on such a quantity. We do not consider its consequences 
on the propagation of the electromagnetic waves involved in ranging. 

Figure [1] depicts the numerically produced range perturbations Ap induced by eq. (JH1) 
for Mercury, Venus, Mars and Saturn over different time spans. The quantitative features 
of the signatures of Figured] are resumed in Table HI It can be noticed that, although their 
temporal signatures are rather distinct, their magnitudes are quite small, amounting to 
about 0.1 — 10 mm. 



3.2. The Moon 

In Table [5] we quote the relevant orbital and physical parameters of the Earth- Moon 
system. Table [6] displays some computed orbital effects of a dipole-type spatial gradient of 
the fine structure constant for the Moon: a geocentric frame with the mean equinox and 
ecliptic at the epoch J2000.0 is assumed in applying eq. (1251) . eq. (127]) and eq. (129]) . 

It may be interesting to notice that the magnitude of the long-term variation of the 



- 14 - 



IHarh-Mercnry range Earh-Vcnns range 




1 2 3 4 5 12 3 4 5 

I (yr) I (yr) 



Fig. 1. — Numerically produced Earth-planet range perturbations Ap, in mm, induced by the 
Stark-like extra-acceleration of eq. ([8]) for Mercury, Venus, Mars and Saturn. We adopted 
the figures of Table [T] for AQ, and eq. ([6]) and eq. §5§ for B and k, respectively. Each 
curve is the difference between the time series for a given Earth-planet range p computed by 
numerically integrating the equations of motion in cartesian coordinates with and without 
eq. OH]). Both the integrations share the same initial conditions retrieved from the WEB 
interface HORIZONS (Author: J. Giorgini. Site Manager: D. K. Yeomans. Webmaster: A. 
B. Chamberlin) by JPL, NASA. They refer to a heliocentric frame with the mean equinox 
and ecliptic at the epoch J2000.0. The total ranges of integration are At = 2 yr for Mercury 
and Venus, and At = 5 yr for Mars and Saturn. 



Table 4: Peak-to-peak maximum amplitude | Ap| max , mean (Ap) and variance o"a p , in mm, of 
the range signals of Figure [1] caused by the Stark-like acceleration of eq. ([8]) for the spatial 
gradient fine structure constant. The total ranges of integration are At = 2 yr for Mercury 
and Venus, and At = 5 yr for Mars and Saturn. 



Planet 


| Ap |max ( mm ) 


(Ap) (mm) 


a Ap (mm) 


Mercury 


23.1 


-0.9 


4.4 


Venus 


14.8 


-0.8 


4.2 


Mars 


60.8 


-2.4 


14.7 


Saturn 


73.7 


-5.1 


16.2 



eccentricity is 10 yr : this implies that the gradient of the fine structur e constant cann ot 
be the cause of the anomalous secular increase of the 



20011: IWilliams fc Dickeyl 120031: Williams fc Boggsll2009[) e 
discussed in literature ( Anderson fc Nieto^OlO : Ioriol 



unar eccentricity (IWilliams et al. 



(9 ± 3) x 10" 12 y^ 1 recently 
201 lh . The predicted perigee 



precession is as large as 3.3 x 10 mas yr ; the present-day accuracy in determining 
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Table 5: Relevant physical and osculating orbital parameters of the Earth- Moon system. 
a is the semi-major axis, e is the eccentricity. The inclination / is referred to the mean 
ecliptic at J2000.0. Q is the longitude of the ascending node, it circulates with a period of 
6798.38 d, and is referred to the mean equinox and ecliptic at J2000.0. u is the argument 
of pericenter: its period is 2191.50 d. G is the Newtonian gravitational constant. The 
masses of the Earth and the Moon are M and m, respectively. AQ = Qg — Q ffi refers to 



the fine structure constant: its value comes from Damour &: Donoghue ( 20 111 ). The orbital 
parameters of the Moon were retrieved from the WEB interface HORIZONS (Author: J. 
Giorgini. Site Manager: D. K. Yeomans. Webmaster: A. B. Chamberlin) by JPL, NASA, 
at the epoch J2000.0. 

a (m) e / (deg) Q (deg) u (deg) GM (m 3 s" 2 ) m/M AQ 

3.81219 x 10 8 0.0647 5.24 123.98 -51.86 3.98600 x 10 14 0.012 -3.2 x 10~ 4 



Table 6: Orbital effects caused by eq. flH]) on the Moon in the case of a dipolar-like gradient 
of the fine structure constant. For the Earth- Moon system the figures of Table [5] were used. 
The variation of the eccentricity and the perigee precession were worked out from eq. (|25|) . 
The position and velocity perturbations were computed according to eq. ( 12T1) and eq. (|2"9~|) . 
respectively. All the effects are averaged over one Moon's orbital revolution. The present- 
day accuracies in determ ining the lunar eccentric ity and perigee secula r variations are of the 



order of 3 x 10 12 yr 1 flWilliams fc Bogga 120091 ) and 10 1 mas yr 1 ([Williams et al.lll996 



Miiller et al.1120071 120081 ) , respectivel y, while the Earth- Moon post-fit range residuals are con- 
sisten tly approaching the mm level (1 Williams et al.l l2004t iMurphy et al.l 120081 ; iBattat et al. 
2009h . 



'(*) 




AR (mm) 


AT (mm) Av R (as) 


Av T (^) 


Av N 


-3 x 10" 14 


3.5 x 10~ 4 


8 x 10- 4 


-6.3 x 10~ 3 0.32 


-0.07 


0.01 



the lunar o rbital precessions is at 0.1 mas yr 1 level ([Williams et al.l 11996c iMuller et al. 



20071 . 120081 ). According to Table EJ the radial and transverse position shifts per orbit are 



of the order of /jm. The radial variation per orbit of the velocity amounts to about 0.1 
mm yr -1 , while the magnitudes of the other two components are one order of magnitude 
smaller. Improvements in both measurement techniques and in dy namical modeling are 
pointing towards post-fit Earth-Moon range residuals at mm level ( jWilliams et al.l 12004 : 



Murphy et al.l l2008t IBattat et al.l 120091 ) . The peculiar temporal patterns of the signals 



investigated may help in separating them from other dynamical standard effects if and 
when the level of accuracy required to determine them will eventually be reached. 



- 16 - 



We numerically integrated the Earth-Moon range over At = 30 yr; its principal 
quantitative features are reported in Table El While the maximum peak-to-peak amplitude 

Table 7: Peak-to-peak maximum amplitude | Ap| max , mean (Ap) and variance <ta p , in mm, of 
the numerically integrated geocentric lunar range signal caused by the Stark-like acceleration 
of eq. (EJ) for the spatial gradient fine structure constant. It was produced as for the planets: 
the total ranges of the integrations areAt = 30 yr. 



| Ap |max ( mm ) 


(Ap) (mm) 


(T Ap (mm) 


4.5 


3 x 10" 4 


0.8 



is as large as a few mm, its mean amounts to 10 4 mm, with a variance of 0.8 mm. They 
are certainly quite small figures, even for a so large time span as the one adopted. 



4. Summary and conclusions 

We considered a dipolar spatial variation of the fine structure constant, for which 
empirical evidence at about 4<r level currently exists, and looked at the impact that the 
resulting Stark-like anomalous acceleration may have on the orbital motion of a test particle 
around a central body. Since the relative two-body Stark-type acceleration depends in a 
certain way on the different composition of the bodies involved, it violates the equivalence 
principle. 

We, first, analytically worked out the long-term, i.e. averaged over one orbital 
revolution, variations of all the six osculating Keplerian orbital elements of the test particle. 
We did not restrict ourselves to any specific orbital configuration; moreover, we made no 
a-priori assumptions on the fixed direction of the gradient of the fine structure constant. 
Thus, our results are quite general. It turned out that, apart from the semi-major axis, 
all the other osculating Keplerian orbital elements experience non-vanishing long-term 
variations which, in real astronomical scenarios, would be modulated primarily by the low 
frequencies of the inclination, the node and the pericenter due to the standard mutual 
N-body interactions. Then, we worked out the changes per orbit of the position and 
velocity vectors of the test particle. Also in this case, the expressions obtained are valid for 
any orbital geometry and gradient direction. Only the long-term shift of the out-of-plane 
component of the position vector vanishes, while the other non-vanishing ones exhibit 
characteristic slow time- varying modulations because of the presence of the inclination, the 
node and the pericenter in their expressions. 

Subsequently, we computed our predicted effects for some Sun-planet pairs of the solar 
system and for the Earth-Moon system. In particular, the perihelion precessions of the 
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inner planets, which exhibit the largest Stark-like accelerations because of their markedly 
different composition with respect to the Sun, are of the order of 10~ 2 — 10 -3 mas cty -1 . We 
numerically computed times series some years long of the Stark-induced perturbation of the 
interplanetary range for some Earth-planet pairs as well; their magnitude is approximately 
0.1 — 10 mm. Concerning the Moon, the long-term precession of its perigee is of the order 
of 10~ 4 mas yr _1 , while the long-term variation of its eccentricity is a decrease at the 10 -14 
yr _1 level. The radial and transverse shifts per orbit of the Moon's position amount to 
some /im, while the radial, transverse and normal variations per orbit of its velocity are 
of the order of 0.1 — 0.01 mm yr" 1 . In addition to that, we also numerically produced an 
Earth-Moon perturbed range time series spanning 30 yr: its main quantitative features are 
at a sub-mm level, with a nominal peak-to-peak amplitude of a few mm and a mean of 0.3 
/im. 

The level of accuracy in empirically determining the corrections to the standard 
Newtonian-Einsteinian perihelion precessions from planetary observations is currently 2 — 3 
orders of magnitude larger than the predicted Stark-like effects, apart from, perhaps, the 
Earth provided that VLBI data are included in the data processing. Recent years have 
seen increasing efforts toward the implementation of the planetary laser ranging technique, 
which should be accurate at cm level. Concerning the Moon, its orbital precessions can be 
presently determined at a 10 _1 mas yr -1 level of accuracy, while the LLR post-fit range 
residuals are approaching the mm level. 

In conclusion, detecting the putative spatial variations of the fine structure constant 
through its orbital effects on the major bodies of the solar system seems unlikely or very 
difficult in the near future. 
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